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A case of explosive extratropical cyclogenesis in the 
North Pacific is analyzed employing quasi-Lagrangian diag- 
nostic techniques in isobaric coordinates. The First GARP 
Global Experiment (FGGE) Level Ill-b data set from the Euro- 
pean Center Fcr Medium-range Weather Forecasting (ECMWF) is 
used in a synoptic investigation and mass and vorticity 
budget diagnostic evaluation of storm development. 

Mid-tropospheric positive vorticity advection (PVA) acts 
in concert with low-level thermal advection to enhance 
surface layer crganizaticn during the initial periods. 
Explosive development occurs under upper-level zonal flow 
when the surface system, developing along the low-level 
baroclinic zone, moves under the front left quadrant of a 
strong jet streak. stability in the lower troposphere 
decreases steadily as an intense two-layer mass circulation 
is established. The divergence term provides the largest 
contributicn to the vorticity tendency during explosive 
development due to substantial surface layer convergence and 
upper-level divergence. 
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I. 



I NTROD UCTION 



Extratropical maritime cylones restrict commercial, 
industrial and military use of the sea. The more intense 
storms which form over very small time periods are espe- 
cially dangerous to navigation due to the effects of the 
associated wind, waves and visibility restrictions. Accu- 
rate and reliable prediction of these atmospheric processes 
has teen a goal for many decades. During the past few dec- 
ades, significant advances in forecasting these synoptic- 
scale atmospheric circulation features have been realized 
through the development of numerical weather prediction 
(NWP) models. 

While current numerical models provide satisfactory 
accuracy in forecasting most synoptic-scale processes, they 
perform poorly in predicting the predominantly cold-season 
explosive genesis of maritime cyclones. Sanders and Gyakum 
(1980) define explosive cyclogenesis as an extratropical 
cyclone which undergoes central pressure falls of at least 1 
mb/h for 24 h. Accurate forecasting of these meteorological 
events is limited by the scarcity of observations over 
oceanic areas, which inhibits the model initialization pro- 
cess, and by the smaller spatial and temporal scales of the 
explosively developing cyclone. 
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Accurate forecasting of these maritime events is of 
great significance to co miner cia 1, industrial and milirary 
organizations which routinely use the oceans in the trans- 
action of business and in the fulfillment of co ounitments. 
Several recently publicized examples, specifically, the less 
of the Russian trawler Metros toy (Mariners Weather Log, 

1979) , and the sinking of the oil rig Ocean Ranger (Lemoyne, 
1982) , illustrate only a minute fraction of the total loss 
of life and operational hazards associated with these 
intense metecrolcgica 1 events. 

Explosive deepening is a character isric of the vast 
majority of the most intense cyclones (Sanders and Gyakum, 

1980) . In connection with associated severe weather and the 
failure of numerical forecast models to accurately predicn 
these events, it is of paramount importance zhat the 
processes which contribute to the formation of these vigor- 
ous storms be examined and fully understood. Several recent 
studies have focussed attention on detailed analyses of 
these processes. Roman (1981) completed a mass budget study 

and Conant (1982) completed a circulation and angular momen- 

\ 

turn budget study of a case of east-coast cyclogenesis. 
In tercomparisens of these results with conclusions derived 
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in this thesis provide valuable insight into the formative 
mechanisms of maritime explosive cyclogenesis and lead to a 
general quantitative understanding of the forecast problem 
associated with these intense systems. 

This thesis is part of a broadly-based research effort 
attempting to understand better the formative processes 
inherent in the genesis of maritime cyclones. Specifically, 
this study concentrates on a clearly defined srrucrured 
approach to the analysis of a case of explosive cyclogenesis 
over the North Pacific Ocean during the period 12-17 January 
1979. The mass and circulation (vorticity) budgets are 
examined using quasi-Lagr angian diagnosrics (QLD) as devel- 
oped by Jchnscn and Downey (1975 a and b) in order to deter- 
mine the role of dynamic forcing in the budget volume. The 
local contribution of vorticity sources and sinks and the 
vertical and horizontal mass and vorticity transports are 
specifically explored. 

Identification of these processes is uniquely interest- 
ing in that the cyclone under study consists of the classic 
frontal wave development from strong low level baroclinic- 
ity. The distinguishing feature is rhe absence of a strong 
migratory short wave aloft, which categorizes this storm as 
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of a cylindrical budget volume centered on and translating 
with the cyclcne surface center. Vertical distributions, 
lateral exchanges and sources and sinics of cyclone proper- 
ties resulting from purely developmental processes are then 
analyzed. Application of QLD techniques in studying cyclone 
development (Mash, 1978) has proven valuable in diagnoses of 
poorly forecast storms, and Roman (1981) and Conant (1982) 
have demonstrated QLD applicability in detailed analyses of 
east-coast cyclogenesis. 

The main objectives of this thesis are; 

• A synoptic investigation of the classical frontal wave 
cyclogenesis in the Kcrth Pacific Ocean; 

• A synoptic assessment of the validity of the First GARP 
Global Experiment (FGGE) data and the quality of the 
FGGE analyses for depicting storm dynamical processes; 

• A diagnostic investigation of storm dynamics through 
the use of QLD budget analyses of mass and vorticity. 

The detailed analyses of these processes is further 
aided by an exceptionally strong supporting daza base. The 
case under study occurred during the first special observa- 
tion period (SOP-1) of the First GARP Global Experiment 
(FGGE). This FGGE data set was assimilated by the European 
Center for Medium-range Weather Forecasting (ECMWF) and 
obtained through the auspices of the Goddard Space Flight 
Center. This data set is special, in that observations from 
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aircraft, ships, buoys and satellites gathered during an 
intense collecrion effort were manually checked and then 
systematically integrated into a dense network of 
observations. Consequently, the supporting data base ana- 
lyzed in this thesis effort is exceptionally broad-based, 
complete, and integrated into a well-diversified 
observational network. 

Quantitative analysis cf FGG E-type data bases gives spe- 
cific insights into development processes inherent in mari- 
time explosive cyclogenesis. In the following chapters, 
examination of the FGGE data (Chapter 2) and an abbreviated 
description cf QLD procedures (Chapter 3) are presented. 
Chapter 4 discusses the cyclone under study through a synop- 
tic overview which highlights significant flow patterns and 
associated upper and lower-level features. Chapters 5 and 6 
present the results of the budget calculations and Chapter 7 
contains specific conclusions from these budget calculations 
and discusses recommendations for further research. 
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II. E NVI RONMENTAL DAT A BA SE DES CRIPTION 
A. FGGE DATA 

Curing the first special observation period (SOP-1) of 
the Global Weather Experiment, the largest concentration of 
meteorological resources ever assembled was deployed to 
observe the atmosphere and oceans. The quality of the data 
generated from this effort marked a milestone in the 
development of numerical weather prediction models and 
synoptic meteorology in general. By any standard, this 
experiment shewed an unparalleled example of scientific 
cooperation and excellent resource management. 

The FGGE Level lll-b data set used in this thesis was 
constructed from observations integrated through a data 
assimilation system at ECMWF, Bengsston et al., (1982) 
describe the system as a three-dimensional, multivariate, 
optimum interpolation scheme using a nonlinear normal-mode 
initialization and an automatic system for data checking. 
The direct field-level observations (Level I-b) were trans- 
formed into the basic meteorological variables (Level Il-b) , 
which constituted the basic meteorological product of the 
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experiment. A 15-level model with horizontal resolution of 
1,875 degrees of latitude and longitude was used to trans- 
form the Level Il-b data to Level Ill-b. At rhe Ill-b 
level, the observations were manually checked and objec- 
tively analyzed into a set of dynamically derived meteorolo- 
gical fields. The flowchart of this data assimilation 
process is depicted in Fig. 1. The FGGS Level Ill-b data 
production effort was completed in June 1981 and the data 
fields were catalogued for 6-h observations for all standard 
levels up to and including 10 mb. A detailed description of 
this assimilation scheme can be found in Bengsston et al., 
(1982) . 

The FGGE data base is currently undergoing validation by 
a variety of groups. Roman (1981) and Conant (1982) demon- 

strated its validity on a case of east-coast cyclogenesis. 
In an independent assessment of the FGGE data set, Halem et 
al., (1982) present a detailed analysis in connection with 
the coverage and accuracy of the remote sensing portion of 
the network covering the oceanic regions. Analyses based on 
these inputs provided the following conclusions; 
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• The FG6E ctserving system has a positive effect 
forecast skill using S1 criteria; 



• Significant corrections to the 300 mb first guess 
fields cccur over oceanic regions with the FGGE data; 



• Nearly all of the forecasts that retain useful skill by 
day five shew a positive impact from the FGGE system; 



• In certain synoptic situations, the FGGE observing sys 
tem is capable of measuring the atmosphere globall 
with sufficient accuracy for current NWP models'to mak 
useful extended- range forecasts beyond one week. 



These and other conclusions as well as a detailed descrip- 



tion of the analysis scheme can be found in Halem et al.. 



(1981) and Baker et al. , (1981) . It is evident from these 



and ether studies that the FGGE Level Ill-b data set pro- 



vides a strong foundation for accomplishing the objectives 



of this thesis. 

A summary of the FGGE data availability by field-type 
and level is depicted in Table I. The advantages of this 
type of augmented data support over the standard data base 
incl ude : 

• A considerably richer data base, particularly over 
data-sparse oceanic areas; 

• The high resolution grid is convenient for deoicting 
the smaller synoptic-scale features associated with 
cyclcgenesis; 



• Basic and derived fields (vertical velocity, humidity, 
etc.) are included in the data set; 
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• The full area of cyclogenesis and srorm translation 
enjoys siotilar data coverage. 

These significant advantages demonstrate the value of the 
FGGE Level Ill-b data set in quantitative analyses of storm 
features, 

B. SUMMARY OF ANALYSES AND PROGNOSES UTILIZED 

Three sets of analyses are utilized in this thesis: 

• National Meteorological Center (NMC) analyses and prog- 
noses; 

• Fleet Numerical Oceanography Center (FNOC) analyses and 
prognoses; 

• ECMWF Level m-b analyses. 

These analyses are used in conjunction with the synoptic 
discussion (Chapter U) and for the verification of the 
forecast problem alluded to earlier. A complete listing of 
all analyses and prognoses utilized in this study is con- 
tained in Table II. A complete description of the FNOC 
products can be found in the U.S. Naval Weather Service 
Numerical Environmental Products Manual (NAVAIR 50-1G-522). 
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III. CUASI-LAGBANGIAN DIAGNOSTIC (^D) 

TEUHTTTU U-E? — 

lo derive generalized budget analysis techniques for 
processes related to cyclone development, changes in the 
budget volume properties due to migration of the system must 
be separated tc isolate true developmental trends. In the 
past, synoptic-scale diagnostic approaches (s.g. Kung and 
Baker, 1975) evaluated storm processes relative to a fixed 
coordinate system. This type of approach is effective for 
selected time periods when storm translation is negligible, 
but suffers markedly when applied to rapidly zranslating 
systems of appreciable intensity. In particular, changes 
may cccur due tc system mcvement, system development, or 
both. The only accurate approach to assess true developmen- 
tal fluctuations associated with cyclone development is tc 
construct a moving budget volume centered on some distin- 
guishing storm feature. In this way, horizontal advection 
associated with wave translation may be isolated. In addi- 
tion, -any vertical redistribution of a budget property may 
be analyzed in terms of vertical transport processes within 
the cyclone budget regime. QLD techniques employed in this 
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thesis provide the means fcr this type of approach. Roman 
(1981) and Conant (1982) have demonstrated rhe validity of 
this type cf approach on a case of east-coast cyclogenesis. 

To study the fluctuations in any storm property, a 
generalized budget equation is constructed in terms of the 
storm budget coordinate system. The equation consists of 
the time rate of change of the desired budget quantity 
(df/dt) and the lateral flux transport, vertical flux 
transport and sources/sinks (see Table III) . The equation 
is applicable to the storm budget volume in spherical coor- 
dinates. The spherical portion of the atmosphere incorpo- 
rated within the budget volume is approximated by a cylinder 
for ease in computations and definition of budget proper- 
ties. This approximation is valid for the relatively thin 
slab of atmosphere (1000-100 mb) selected for budget 
anal yses . 

The formulation of the storm budget coordinate system 
and details cf the budget equation notation are discussed by 
Hash (1978). The storm budget volume is defined in terms of 
10 mandatory pressure levels in the vertical, and in the 
horizontal by outwardly directed radii originating from a 
central point. For example, the budget property at a point 
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may te analyzed by describing the level/layer (1000/925 mb), 
the budget radius in degrees of latitude from the storm 
center (6) , and the quadrant desired (grid points 1-36) . 
The mass and vorticity budget equations employed in this 
analysis are depicted in Tables IV and V respectively. 

In the mass budget (Table IV) , the general equation rep- 
resenting the change in mass with time is composed of con- 
rributions from lateral and vertical transports. The 
lateral transport term represents the ner horizontal conver- 
gence of mass flux within the budget volume. The vertical 
transport is calculated through the use of derived vertical 
velocity fields and is necessarily a reflection of ner 
lateral flux into the volume. Application of boundary con- 
ditions that vertical motion is zero at the top and bottom 
surfaces of the budget volume simplifies rhe equation, as 
shown in the middle of Table IV. 

Several useful relationships can be derived from this 
simplified equation and the boundary conditions. Contribu- 
tions to the mass time tendency in iso baric coordinates 
arises chiefly from pressure fluctuations in the surface 
layer, as the other layers in the budget are fixed pressure 
surfaces. This term is compuxed easily and with 
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In the vorticity budget formulation, the contribuxior.s 
to vorticity changes occur from lateral and vertical fluxes, 
sources and sinks (Table V, section A) . The flux 
(transport) terms can be expanded into divergence and advec- 
tion contributions by the use of vector identities (Table V, 
Section B) , which facilitates further insight into the phys- 
ical processes they represent. The source/sink terms arise 
from net generation/dissipation of vorticity within the 
budget volume, and in iscbaric coordinates are comprised of 
divergence, tilting and friction components. Notice that 
three divergence terms arise from this expansion of the 
general vorticity budget equation. The horizontal transport 
component is identical to the source divergence term, and 
the horizontal and vertical transport components are equal 
and opposite. Thus, the net contribution of divergence in 
the budget is manifest in the source term. A specific 
breakdown of the partitioned vorticity budget equation is 
presented in Section B of Table V. 

To isolate the contributions of both the rotational and 
irrotational wind components to the transport of a desired 
budget quantity, mean and eddy modes (analogous to 
divergence and advection terms in the vorticity equation) , 
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are derived from the transport terms (Table VI) . The pur- 
pose of this technique is to isolate the changes in the mean 
budget properties from these occurring on the eddy scale. 
Thus, a comprehensive analysis of the lateral and vertical 
flux terms can be made by separating them into mean and eddy 
mode contributions and comparing them with the divergence 
and advection terms for analysis of the physical processes 
they represent. The relation of the mean/eddy mode and 
divergence/advection gartitions to the total transport is 
also illustrated in Section B of Table V. 

The budget analysis begins with the choice of the par-<- 
ticular property to be analyzed. The next step in the 
budget calculations consists of interpolating the u and v 
compenents of the ?GGE data sat onto the cylindrical budget 
volume grid. Next, the interpolated winds are converted to 
tangential and radial compenents relative to the outer vol- 
ume boundary and corrected for storm motion. In this form, 
the data are ready to be used in the budget equations 
discussed above. Derived fields (moisture, vertical veloc- 
ity, kinetic energy, etc.), not directly used in budget com- 
putations can be used to supplement insight into storm 
dynamics. Interesting comparisons between FGGE derived 
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vertical velocity and QLD calculated vertical velocity, for 



example, serve to validate the QLD technique approach. 



The QLD approach provides insight into the following key 



properties : 



• The magnitude and significance of the horizonral mass 
transport structure; 



• The magnitude and significance of the horizontal vor- 
ticity transport structure; 



• The vertical transport structure of mass and vorticity 
derived from vertical motion estimates; 



• The magnitude and nature of source/sink contributions 
to budget analyses. 



The next chapter discusses the development sequence of the 



storm through a synoptic discussion urilizing synoptic 



fields, analyses and forecasts, and satellite imagery. 
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IV. 



SYNOPTIC OVERVIEW 



A. GENERAL 

Oceanic cyclogenesis mechanisms can in part be portrayed 
by graphical representation of pertinent meteorological and 
oceanic fields. This chapter explores the synoptic condi- 
tions associated with a specific case of North Pacific 
explosive cyclcgenesis during 12-17 January 1979. An 
attempt is made to gain insight into formative processes and 
possible factors related to the serious shortcomings in the 
model repesentaticn cf the dynamics associated with these 
vigorous events. Use is made of FGGE data analyses. Defense 
Military Satellite Program (DMSP) imagery and FNOC sea 
surface temperature (SST) analyses. Comparisons are mads 
between FGGE , FNOC and NMC sea level pressure analyses and 
prognoses to highlight differences in forecast deepening 
rates and analyzed and forecast storm positions. Explosive 
deepening criteria were adopred from Sanders and Gyakum 
(1980), who addressed the NMC 7-layer PE model performance 
associated with "explosive" prediction. Predicted deepening 
rates for the cyclones they studied were consistently less 
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than 1/3 of the observed rates, and North Pacific explosive 
deepening east of the dateline was virtually missed by the 
model. An assessment of the FNOC 5-layer PE model perform- 
ance (Gyakum, 1981) demonstrated similar results. 
Specifically, the model captured only 1/3 of the observed 
central pressure falls for all time periods. In addition, 
the mcst accurate of all the forecasts were shown to be 
valid between 24-36 h. In comparing the FGGE analyses with 
the FNOC and NMC analyses, it should be noted that the FNOC 
analyses and prognoses are derived from the 5-layer PE 
coarse mesh (381 km hcrizcntal resolution) model. In con- 
trast, the NMC analyses are strictly hand-analyzed Northern 
Hemisphere strips, while the prognoses are derived from the 
7-r layer PE fine mesh (190.5 km horizontal resolution) model. 

B. MODEL FREEICIION PERFCEMANCE 

The storm track (at 12-h increments) is presented in 
Fig. 2a. The cyclone develops southeast of Japan in the 
vicinity 30N, 135E and transits toward a final position off 

the west coast of Alaska. Notice the significant transla- 
tion increase after 0000 GMT 13 January. The 5-day evolu- 
tion (0000 GMT 12 January 1979 to 0000 GMT 17 January 1979) 
in central sea-level pressure, as analyzed from the FGGE 
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and t.his caused a misrepresen- 



closest to the storm center, 
ration of the surface central pressure. Thus, the more 
accurate representation of the deepening phase in the FGGS 
analysis is directly attributed to a superior data base and 
not just to the benefit of a "better" first guess. 

A comparison between the 36-h sea-level pressure 

forecasts is shown in Fig. 3a. All curves show the general 

deepening trends, but the FNOC and NMC forecasts 

consistently lag the FGGE analyzed pressure, and are highly 
variable in predicting the magnitude of the central pressure 
falls. From these displays, it is evident that both models 
fail significantly in the timely and accurate prediction of 
storm intensification, especially during the explosive 
stage. This example illustrates further the forecast prob- 
lems discussed by Sanders and Gyakum (1980) alluded to 

ea rlier. 

The location error in 36-h storm positions for NMC and 
FNOC predictions, relative to the FGGE analyzed position, is 
presented in Fig. 3b. It is clearly seen that significant 
errors exist in both forecasts at various times curing the 
deepening phase. Such significant errors in the analyzed 
and forecast features of this system reflect an inability of 
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current operational models to accurately forecast these 
events. 

C. SYNOPTIC DISCUSSION 

The synoptic analysis approach in this section consists 
of an examination of the basic meteorological fields in 
three displays of two fields each. Display A contains 250 
mb winds and isotachs, B presents 500 mb heights and abso- 
lute vorticity, and C shows 1000-500 mb thickness and sea- 
level pressure patterns. This approach generally allows 
only subjective insights into features of interest during 
the storm life cycle. The quantitative analysis of the 
dynamical processes is presented in the following chapter. 
The intent of this section is to trace the development of 
the system through a syncptic discussion of major features 
contributing to overall storm development. 

Examination of atmospheric conditions at 0000 g:iT 12 
January 1979 in the area southeast of Japan shows marked 
cyclcgenetic potential. The surface picture (Fig. 4c) is 
dominated by a broad trough (vicinity 30N, 140E) extending 

northeast-southwest through the center of a checkerboard 
pattern of fairly strong low and high centers. The signifi- 
cant features include the Siberian High near 55N, 100S (jusr 
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off the western edge of the chart) which is very strong and 
stationary. Ihe parent lew in the vicinity of 53N, 170E is 
now quasi-stationary and filling (950 mb) . A large maritime 
high (1030 mb) dominates the north central Pacific (vicinity 
30N, 178W) , and the northern fringe of Typhoon Alice is jusn 
visible near 20N, 1U0E) . The combination of flow from rhese 
systems serves to funnel cold, polar continental air south- 
ward towards northward flowing warm, moist tropical air, 
which enhances the frontogenes is potential. A general 
destabilization of the polar air occurs norrh of the area of 
strong low level baroclinicity shown by the concentration of 
thickness contours just east of Japan. 

In the mid-troposphere (Fig. 4b) , two distinctive fea- 
tures are evident. The flow over the western Pacific is 
dominated by an extensive region of low heights associated 
with a quasi-stationary long wave trough, which is a portion 
of a four-wave planetary pattern. This significant feature 
accounts fer the steady progression of mobile short waves 
and associated surface features through the region during 
the preceding week, in addition to supporting the steady 
flow of cold air off the Siberian continent. The second 
major feature of interest is the predominantly zonal flew 
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over ‘he entire incipient region with areas of weak and 
small absolute vcrticity upstream. 

In the upper troposphere (Fig. 4a), the presence of a 
broad jet ( > 50 m/sec) aloft is a further manifestation of 
the baroclinicity of the lower troposphere. As yet, signif- 
icant jet streak ( > 75 m/sec) propagation into the region 
remains well upstream, and therefore is not a factor in the 
low-level development. 

Examination of the visual (Fig. 5) and infrared (Fig. 6) 
satellite imagery valid at 0041 GMT 12 January 1979 shows a 
couplet of similar low level cloud masses (labelled A and B) 
southeast of Japan. Lack of significant vertical 

development is evident on the infrared imagery, in contrast 
to the area of enhanced convection (labelled C) associated 
with an area of positive vcrticity advection (?VA) upstream 
over Korea (Fig. 4b) . Cloud mass B is the precursor of a 
weak cyclone which moves to the east and rapidly dissipates. 
Cloud mass A is the precursor of the cyclone analyzed in 
this study. The evolving low-amplitude "s-shape” configura- 
tion with a distinct edge on the western flank is indicative 
of increasing thermal and low level absolute vorticity gra- 
dients with time (Weldon, 1977) . Although fairly small, the 
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cloud mass is beginning tc take on the familiar "leaf -shape" 
which Weldon (1977) characterizes as rypical of incipient 
oceanic cyclogenesis. 

During the next 24 h ending at 0000 GMT 13 January 1979, 
surface development continues (vicinity 35N, 150E) into a 

closed center (995 mb). The incipient sysnem (Fig- 7c) 
appears as a broad area of low pressure, and the development 
of a weak thermal ridge near the surface center indicates 
frontogenesis continues during the period. Ships in the 
vicinity reported seas on the order of 15 ft, with periods 
ranging between 3-5 s, which indicates a duration limited 
wave spectrum. Of interest is the continued strengthening 
of the Siberian high (now 1080 mb) which steadily pumps cold 
air southward, in association with the backside circulation 
of the stationary cold low to the east. Dasrabilization of 
the surface layer continues as cold air streams southward 
over the warm ocean surface, and a polar low has become 
organized ever the southeast tip of Japan due west of the 
incipient sysxem (indicated on detailed NMC analyses, not 
shown here) . 

In the mid-troposphere (Fig. 7b) , the 500 mb low 
anchored off the Kamchatka Peninsula has deepened. A short 
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wave with at associated broad absolute vorticity area has 
formed over Japan, approximately 800 nm upstream of the 
incipient system. A smaller, less distinct vorticity area 
has formed directly over the surface system. 

Increased vorticity in the mid-troposphere can also be 
linked to the approach of a jet streak from the west. At 
this stage, however, the position of the jet streak indi- 
cates only minimal influence on the low-level circulation, 
as the jet core is directly superposed in the vertical, and 
no clear pattern of cyclonic shear or diffluence can be seen 
at 250 mb (Fig. 7a) . 

The increased vertical development over the period can 
be seen in the DMSP satellite imagery (Figs. 8 and 9) valid 
at 23U 1 GMT 12 January 1979. The developing system is 
labeled A and is centered approximately 600 nm southeast of 
Japan. Cold plumes extend for virtually the entire breadth 
in the infrared image (Fig. 9) , with cloud tops becoming 
progressively warmer toward the tapered end to the south- 
west. Enhanced convection and continued development of the 
leaf structure is shown in the visual image (Fig. 8) . The 
region under the leaf appears to contain marked baroclinic- 
ity in a zone in advance of the propagating jet streak 
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alcft, which is similar tc the conditions observed by Weldon 
(1 977) in a case study. 

Curing the following 24-h period ending 0000 14 Jan- 

uary 1979, continued development is evident in nhe low and 
middle troposphere. Central surface pressure (Fig. 10c) 
falls steadily tc 977 mb, which is a decrease of 18 mb. If 
one adopts the explosive criteria of Sanders and Gyakum 
(1980) for a latitude of 40 degrees, rhe srorm has entered 
the explosive stage during this period. The surface low 
moves rapidly toward the northeast, to a position close to 
the maximum SSI gradient (Fig. 11). Close proximity zo this 
gradient fulfills another observation of explosive 
development described hy Sanders and Gyakum (1980). 
Analysis of the thickness contours (Fig. 10c) indicates a 
pattern of diffluence 10 degrees of longitude downstream 
from the surface center. The diffluence at 10 degrees 
upstream is ccnsiderably greater. Bjerknes (1954) shewed 
that these conditions favor formation of a diffluent trough 
in the upper levels within the cyclonic shear zone north of 
the jet streak axis. The NW-SE tilt of the 500 mb trough 
evident by 1200 GMT 14 January (Fig. 16b) , combined with the 
favorable superposition cf the jet streak support the 
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diffluent trough formation argument. In addition, a combi- 
nation of the super posit ion of the upper level divergence 
zone over the surface center and the approach of the upper- 
level thermal trough in close proximity to the 500 mb trough 
promotes self-amplification through the mid-troposphere. 
Mid-level vorticity increases due to the movement of concen- 
trated cyclonic shear vorticity into the budget volume as 
the jet enters the diffluent trough. As a result, increased 
vertical development occurred from vorticity increases as a 
result of jet streak interaction and self-amplification. 

Indeed, the first evidence of mid-tropospheric 

development is visible in Fig. 10b. The vorticity pattern 

has an expanded area of significant vorticity associated 
with the mid-tropospheric circulation. A region of strong 
vorticity also exists upstream northeast of Japan, but -hat 
center appears tc be associated with the polar low forming 
to the west (Fig. 10c), as indicated earlier. The interest- 
ing feature tc notice is the absence of the short wave that 
was present ever Japan 24 h earlier, which would seem to 
indicate that further explosive development from the short 
wave PVA aloft is unlikely. 
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kt 250 mb (Fig. 10a) , however, the arrival of the jet 
streak (maximum velocity 85 m/sec) in favorable superposi- 
tion over the surface system acts to promote the lover-level 
circulation. Diffluence in the region of the jet stream 
split and cyclonic wind shear to the left of the jet axis 
combine to induce positive circulation increases and strong 
vertical motion through the mid-troposphere. Increased 
vertical motion enhances vertical redistributions of all 
properties, including an increase in the convergence in the 
surface layer, and hence, explosive development continues. 

Inspection of satellite imagery (Figs. 12 and 13 valid 
1405 GMT 13 January 1979 and Figs. 14 and 15 valid at 0023 
GMT 14 January 1979) shows a discernable comma-shape 
formation. In Figs. 12 and 13, the southern boundary of the 
jet axis is delineated by subsidence along the cirrus band 
paralleling the low-level front. The jet crosses just sourh 
of the vortex at A. The very ccld tops evident near the 
vortex center and along the frontal band indicate strong 
convection throughout the region. An interesting feature in 
Figs. 14 and 15 is the relatively low level clouds (labeled 
A) that appear as fingers extending from the vortex center 
and look similar to feeder bands. When coupled with the 
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extensive cirrus shield aloft, the cloud pattern presents a 
hurricane-like appearance, Twc significant changes have 
occurred during the past 10 h. The area of deep convection 
has expanded markedly both near the vortex and along the 
frontal band, which indicates continued strengthening. 
Second, the storm has organized to the point that clear 
slots have formed in the cold air behind the front, which 
indicates the penetration of cold and subsiding air toward 
the vortex. 

Turing the next 12-h period ending 1200 GMT 14 January 
1979, the surface pressure falls an additional 11 mb to 965 
mb {Fig. 16c). Notice the shift of the reference grid 10 
degrees to the north and 30 degrees to the east in this and 
subsequent figures. The increased amplitude of the 

thickness ridge implies increased circulation in the lower 
troposphere, most probably as a result of warm advection 
and extreme latent heating from deep convection. Sensible 
heat input presumably weakens during the period as the storm 
center moves north of the maximum SST gradient (shewn 
earlier in Fig. 11). 

At 500 mb (Fig. 16b) , self- amplif icat io n continues with 
the development of the upper-level trough and associated 
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vorticity maximum. The apparent lack of appreciable vortic- 
ity advection at this level suggests that continued intensi- 
fication of the system must be manifest in either the lower 
or upper levels. The mcst likely factors which serve to 
maintain the deepening phase consist of a combination of 
processes: jet stream interaction an 250 mb; and the 

extreme surface convergence generating pockets of sustained 
deep convection. 

ftn increase in cyclonic vorticity appears at the 250 mb 
level during this time frame (Fig. 16a). Cyclonic curvature 
and shear to the left of the jet axis continues to promote 
the circulation discussed above. The wind speeds in the jet 
have diminished slightly in intensity due to the decreased 
baroclinicity through the lower troposphere. 

Strong development of the system during the previous 12 
h is indicated in Figs. 17 and 18 (valid at 1205 GMT 14 Jan- 
uary 1979). Vertical velocities within the storm volume 
increase during the period, as deep convection becomes 
intense and widespread surrounding the vortex and along the 
entire frontal band. The growth of convective activity 
indicates that development continues upward through the 
troposphere during the period, as the jet crosses the 
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frontal band at a more northerly point. In addition, the 
"hook" or tip of the cloud band in the northwest sector 
develops markedly. In typical cases (Weldon, 1977) , these 

conditions indicate the onset of the mature stage. 

Interestingly, the next 6-h period ending 1800 GtlT 14 
January 1979 marks the most explosive period of the storm's 
development cycle. The surface pressure (Fig. 19c) falls 18 
mb (averaged between FGGE and NMC analyses). Ship 

observations at 1200 and 1800 GMT were limited in both 
analyses. It is possible that seme of the extreme deepening 
observed during this 6-h period actually occurred before 
1200 GMT, but there were inadequate ship observations to 
detect it. Thus, the magnitude of the surface pressure fall 
may be less than mentioned above or presented in Fig. 3b. 
In any event, the nature of the explosive deepening has 
caused the thickness ridge to double in amplitude along an 
east-west orientation, indicating that cold air has begun to 
wrap around the vortex at the surface. In the mid- 

troposphere (Fig. 19b), the positive vorticity center elon- 
gates along the flow and intensifies slightly. The 

elongation of the vorticity center and the continued digging 
of the mid-level trough indicate the system should weaken in 
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00 ?.^ inuanc€ cf strong vertical motion induced by 



deep con- 



vection alcng the frontal band, and to the north and east of 
the vortex center. The penetration of cold air into the 
vortex is visible and signals the start of the occlusion 
phase . 

At 0000 GMT 15 January 1979, the surface circulation is 
dominated by the surface lew over the Aleutian Islands (Fig. 
22c) with a central pressure at 947 mb. The explosive 
period has ended with the onset of the occlusion process. 
The thickness gradient has diminished markedly, and a polar 
low is organized in the cold air to the west, while a new 
disturbance forms alcng the baroclinic zone to the south- 
west. At 500 mb (Fig. 22b) , a closed circulation has 
formed, which indicates cold air influx through most of the 
troposphere. A broad positive vcrticity pattern exists, but 
negative vorticity advection into the region is inducing 
subsidence and filling. The system is vertically stacked, 
which also indicates the fully mature stage. At 250 mb 
(Fig. 22a), development cf the trough is evident from the 
cyclonic rotation of the wind vectors. The jet has dimin- 
ished in intensity due to the reduced low-level baroclinic- 
ity and the migration of the surface system north of the jet 
ax is . 
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During the next 24 hours ending 0000 GMT 16 January 
1979, the occlusion process continues. The surface systei 
(Fig. 23c) has stalled over the west Alaskan coast and 

filled to 955 mb. Ship reports in the vicinity show 35 ft 
seas with 30 ft swell and indicate that the storm is still a 
potent force. Other signs indicate, however, that the storm 
is beginning to dissipate. The decrease in the thickness 
gradient indicates that the lower tropospheric front is 
occluded. The 500 mb absolute vorticity (Fig. 23b) decrease 
indicates a large reduction in relative vorticity from neg- 
ative vorticity advection and the discontinuance of rhe jet 
stream influence (Fig. 23a). Extreme height falls are 
reflected in the 250 mb flow as the wind field becomes 
closed, indicating an upward cyclonic development through 
this level. 

Inspection of the satellite imagery for this period 
(Figs. 24 and 25) valid at 1005 GMT 15 January 1979 shows a 
fully mature occluded system. The penetration of cold air 
to the vortex center is evident along with a marked decrease 
in deep convective activity near the vortex and along the 
frontal band, except in the northern fringe where the strong 
cyclonic flow is still producing some low-level convergence 



46 



and associated deep convection. The developing low zo the 
south is just visible in the imagery. 

By 0000 GMT 17 January 1979, the surface system has 
filled to 962 mb and the thickness gradient is small (Fig. 
26c) . At 500 mb (Fig. 26b) , the low has combined with a 
stationary feature to the west to produce a vast area of low 
heights. Positive vorticity has decreased as the vertical 
gradients have deteriorated markedly while the system slowly 
dies. At 250 mb (Fig. 26a), the increased closed 
circulation indicates vertical stacking through the whole 
troposphere. These features suggest that the storm has 
reached the dissipating phase, which marks the termination 
of the syncptic portion of the analysis. 

It is interesting to notice that during the period of 
lowest central pressure (0000 GMT 15 January 1979), ship 
reports in the region indicate 40 ft waves, 20 ft swell and 
38 kt winds. The dramatic deepening of this and similar 
storms is net without its consequences. These storms pose 
significant threats to ships and land facilities within 
their paths. The intent of the next chapter is to pursue 
those processes qualitatively described in this chapter by 
mass and vcrticity budget analyses. The budget framework 
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V. 



MASS BODGET A NALYS IS 

A. GENERAL 

It is cl€ar from the synoptic discussion (Chapter 4) 
that the key factors necessary for cyclone development are 
inherently related to dynamical processes occurring in the 
lower tropospheric regime. One such vital process is the 
dest abili2aticn of the surface layer. Another includes the 
horizontal and vertical mass circulation structure, as pres- 
sure falls and resultant storm intensity are directly depen- 
dent upon the nature and magnitude of the mass flux into and 
out cf the budget volume. 

The mass budget portion of this chapter includes a guan- 
titative approach to determine the destabilization and hori- 
zontal and vertical mass transport within the budget volume. 
The stability analysis is based on area-averaged potential 
temperatures fcr the layer. For example, the layer labelled 
775 mb represents the 850-700 mb layer. The lateral mass 
transport is also calculated for the layer and is averaged 
over time periods. Vertioal mass transport is inferred from 
the vertical velocity (omega) fields. Time periods refer to 
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the 12-h tine increment between synoptic times. For exam- 
ple, period 1206 refers to 0000-1200 GMT, Radial distances 
of 6 and 10 degrees of latitude specify the horizontal 

budget area and correspond to 360 and 500 n mi radial dis- 
tances from storm center. These particular radii were cho- 
sen after inspection of radial sections depicting the 

storm's major circulation features. The inner radius (6) 
incorporates the zone of maximum cyclone circulation associ- 
ated with the strong two-layer i nf low/cutf low regime (subse- 
quently addressed). The outer radius (10), in contrast, 
shows the weaker storm circulation near its periphery. 

B. POTENTIAL TEMPERATORE AND STATIC STABILITY ANALYSIS 

The steady flow of polar continental air from the Sibe- 
rian mainland over the warm waters of the western Pacific 
was described in Chapter 4. Destabilization of the boundary 
layer in this manner is a critical factor for oceanic cyclo- 
genesis, as demcnstrated by Roman (198 1) and Sandgathe 
(1981). The amount of destabilization is contingent upon 
the latent heat release in the layer and the sensible heat 
flux from the ocean. These quantities are dependent to a 
large degree on the air-sea temperature difference, and the 
magnitude of the lapse rate through the lower troposphere. 
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Roman (1981) found stability decreases of nearly 2 K / 100 
mb fully 12 h before the cnset of maximum deepening using a 
budget volume analysis. Numerical model simulations by 
Sandgathe (1981) showed a similar trend, with magnitudes 
roughly half that observed by Roman (1981) for a western 
Atlantic storm. An analysis of the destabilization of the 
lower troposphere in this case will be described using the 
potential temperature vertical-time sections. 

Inspection of the time sections of area-averaged 
potential temperature for radii 6 and 10 (Fig. 27a and b) 
shows similar interesting features. Notice that potential 
temperature values in this figure are presented for 1 2-h 
synoptic times. Lower tropospheric (1000-750 mb) and upper 
tropospheric (450-175 mb) potential temperature changes 
imply warming during the early time periods. In the mid- 
troposphere (750-450 mb) , potential temperature changes lit- 
tle in the early time periods. Of interest is the apparent 
lack of a well-defined trcpopause (indicated by a dense 
vertical packing of isentrcpes) during the early time peri- 
ods. Possible factors accounting for this observation 
include the impact of jet streak propagation (described in 
Chapter 4) into the budget volume, or the characteristic 
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upward slope of the tropopause toward the tropics, or a com- 
bination of both. At later time periods, the entire budget 
volume is seen to cool, with the most significant decrease 
occurring below 350 mb, as the cyclone takes on familiar 
ex tratrop ical cold-core characteristics. The tropopause 
becomes well established by 0000 GMT 15 January, which 
coincides with the movement of the budget volume cut of the 
region of the fcaroclinic zone and into the colder polar air. 
Several substantial potential temperature fluctuations exist 
during the 1800 GMT times (not shown) and are addressed 
later. 

To convert potential temperature changes to static sta- 
bility features, a procedure similar to that employed by 
Sandgathe (1981) will be used. He defines the static sta- 
bility of the layer as directly proportional to the specific 
volume of the layer and to the change in potential tempera- 
ture between pressure surfaces (d©/dP) , and inversely pro- 
portional to the average potential temperature of the layer. 
For periods of interest here, the static stability of the 
1000-500 mb layer will be explored. The specific volume of 
the layer is assumed constant with time (variations are less 
than 3%) . The potential temperature fields (Fig. 27) imply 
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rather uniforni average temperature structure for the lower 
troposphere through time. Thus, the static stability varia- 
tions are predominantly a function of the change in the 
lapse rate with time. 

The static stability trends for inner and outer radii 
are shown in Fig. 28a. Destabilization of the layer occurs 
for bcth radii until 1200 GMT 14 January. The increase in 
stability after 1200 GMT 14 January is a reflection of cold 
advection in the lower troposphere. The mariced decrease in 
stability during the early periods (Fig. 28a) is associated 
with an increase in the surface layer potential temperature. 
Thus, surface instability promotes the generation of organ- 
ized convection in the stcrm region, and serves to foster 
the cyclonic circulation tendency. A similar trend was 
observed by Sandgathe (1981) in his numerical simulations 
(Fig. 28b). Direct comparison of observed and simulated 
stability magnitudes is not feasible due zo -he significanr 
differences in the horizontal scales over which the data was 
averaged. It can be seen, however, that rhe trends in the 
observational data match the simulated data very closely, 
and show the characteristics of one of the main ingredients 
that is favorable for oceanic cy do genesis . 
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0ns interesting feature of Fig. 27a deserves special 
scrutiny. The dramatic fluctuations in potential tempera- 
ture (8 K) between 1200 GMT 14 January and 0000 GtIT 15 Janu- 
ary between 800 and 500 mb is remarkable for such a small 
time interval. In general, the atmosphere requires much 
longer periods fcr such large fluctuations. Inspection of 
850-500 mb potential temperature fields for both radii for 
times immediately preceding and subsequent to 1200 GMT shows 
significant warming evident in the 1800 GMT field. In addi- 
tion, potential temperature biases of smaller magnitudes 
were evident at the 1800 and 0600 GMT times on other days 
(Fig. 27a and b) . The nature of these (1800 and 0600 GMT) 
potential temperature fluctuations suggests a problem with 
the FGGE data assimilation schemes, and its impact on this 
study is specifically addressed in Appendix A. 



C. HORIZCNTAI MASS FLUX 
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Inspecticn of Fig. 29a for radius six shows a rwo-layar 
regime consisting of strong low-level mass influx associared 
with surface layer convergence below a layer of mass outflow 
associated with middle and upper-level divergence. The 
strength of the mass flow is seen to increase in both 
regimes as the storm develops, with a distinct inflow maxi- 
mum near 925 mb and an outflow maximum near 350 mb. The 
inflow maximum (3.27 X 10 * * 13) at 1318 GMT is roughly 

twice the outflow maximum (1.56 X 10 * * 13). Although the 
low-flevel inflow is intense, it is restricted to a rela- 
tively narrow layer. The deeper, though weaker, upper-level 
outflow causes a net vertically integrated mass loss within 
the column. The net horizontal mass outflow in the vertical 
columns (Fig. 29a) is reflected in the surface pressure 
falls (Fig. 2t) during the developing periods. 

The temporal variation in the level of non-divergence 
(LND) and the associated inflow maxima (Fig. 29a) between 
1306 and 1406 GMT indicates vertical growth of the 

convergent layer as intense cyclogenesis occurs. By 1306 

GMT, when closed isobars become evident at the surface, the 
convergent layer begins tc grow as the surface circulation 
is established. The growth of the convergent layer during 
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cyclcgenesis is consistent with results found by Honan 
(1981). The radius 6 LND reaches maximum height (at 1406 
GMT) soon after the occurrence of the maximum surface layer 
transport and just before the period of lowest surface pres- 
sure (1500 GMT), and then exhibits a gradual reduction with 
time as the stcrm begins to dissipate. 

Comparable features exist within radius 10 (Fig. 29b) 
but interesting differences exist. The mass buildup in the 
outer volume (1318 GMT) shows a greater relative increase, 
which indicates differential transport between the two 
radii, since the inner vclume increases are incorporated in 
the cuter volume tendency. The differential mass transport 
can be explained by the earlier influence of the migratory 
jet streak (discussed in Chapter 4) as it propagates first 
into the outer budget volume. The decrease in mass 
transport after 1318 GMT reflects the decrease in jet streak 
influence as it propagates through the inner volume and dis- 
sipates (see Chapter 4). These patterns substantiate the 
intensity cf associated '•weather'* observed by stations hun- 
dreds of n mi from the stcrm center and clearly illustrate 
the immense scale of these meteorological events. It should 
also be noticed that the increase in low-level circulation 
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during the last few time periods at radius 10 (Fig. 29b) is 
attributed tc the incorporation of the developing 
circulation to the south (see Chapter 4) . 

Several interesting features of the development stage of 
this cyclone can be observed by inspection of the radial- 
vertical sections of the corrected lateral mass flux. A 
sequence cf fcur times was chosen (0000 and 1200 GMT 14 Jan- 
uary and 0000 and 1200 GMl 15 January) illustrated in Figs. 
30 and 31 respectively. 

At the outer radii in Fig. 30a, the mass transport is a 
maximum near radius 11. The LND is relatively constant in 
height through all radii. A broad region of strong upper 
level divergence is also evident near the outer radius. 
Notice the significant change twelve hours later (Fig. 
30b) . An additional layer of surface convergence near 775 
mb approaches the storm center from the outer radii. This 
feature suggests that a major significant change in the 
circulation has occurred during the past twelve hours. Nor- 
ma] cyclone development might explain a portion of the 
increase, but the explosive development indicated during 
this time from the rapid pressure falls suggests some 
additional forcing present in the budget volume. 
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In the subsequent 12-h period (Fig. 31a) , the surface 
transport maxiniuni decreases and moves vrell inside radius 6. 
Divergence aloft is well established and the positive 
transport inside radius 3 between 400-125 mb suggests net 
subsidence and adiabatic warming near the cyclone core, 
which arrests at least to some small degree the cold-ccre 
formation process. Transport values markedly decrease 12 h 
later (Fig. 31b) and suggest cyclone dissipation. 

These circulation features may be compared with major 
characteristics of tropical storm intensification (e.g. Hol- 
land and Merrill, 1983). Explosive development in both 
cases appears to be characterized by large circulation and 
mass transport increases at the outer radii, with subsequent 
inner radial development after organized transporr has moved 
near rhe inner region with deep convection. The inner 
radial development is also largely dependent on the 

formation of a broad and strong outflow regime in the upper 
troposphere. As these significant features form, organized 
convection dominates the inner radii and triggers the rapid 
pressure falls due to large latent heat inputs into the 
lower troposphere. Large pressure falls set the stage for 
increased winds, circulation etc.. The storm rapidly 
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develops with warming near the central core until the 
surface system becomes detached from the upper-level venti- 
lation and divergence source. The significant observed dif- 
ference in the formation processes lies in the expanded 
radial scale of the extratropical system compared to the 
mesoscale structure of the tropical system. Otherwise, 
upper-level ventilation features, inner radius organized 
deep convection, weak stability of the surface layer and 
large outer region transpcrts in the early development peri- 
ods suggest reasonable similarities. The conclusions cited 
above are only suggestions from the available FGGE fields 
and will not be explored in detail, as these aspects are not 
considered a primary objective of this study. 

D. VERTICAL MASS TRANSPORT 

The storm vertical mass transport analysis consists cf 
vertical time sections of cmega fields (Figs. 32, 33 and 34) 
for radius 6 and 10. Three distinct omega fields are calcu- 
lated by slightly different methods and rhe differences are 
examined. The omega fields discussed below represent the 
vertical mass transport and necessarily reflect periods of 
intense hcrizcntal mass transporr through continuity (see 
Chapter 3) . 
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The differences between the three omega fields are shown 
clearly in Figs. 32, 33 and 34. The budget calculated 

(QLD) omega fields (Fig. 32) are area-averaged values 
derived from the lateral transport through mass continuity. 
The FGGE Level Ill-b vertical motion estimates (Fig. 33) are 
derived as point estimates without diabatic internal adjust- 
ments. Omega fields calculated at grid points using the 
kinematic method, which implicitly includes small-scale 
diabatic effects, are shown in (Fig. 34). 

The difference in magnitude and level of the vertical 
motion maxima between the three vertical motion fields is 
interesting. The QLD (Fig. 32) and FGGE (Fig. 33) vertical 
velocity maxima match fairly closely in time of occurrence 
and magnitude, but differ considerably in level. Compensa- 
tion for the residual imbalances in the budget calculations 
(discussed in Chapter 3) are applied layer by layer (highest 
layers are weighted more) in the QLD omega field derivation. 
Thus, the QLD omega magnitudes are comparable but the level 
of the maxima is biased toward the higher layers (Fig. 32). 

The kinematic vertical velocities differ both in magni- 
tude and level. This is not surprising considering the val- 
ues include small-scale diabatic effects. Paegle (1983) 
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demonstrated similar results and found that the FGGE derived 



fields were roughly 1/3 the magnitude of the kinematically 
derived fields, which is consistent with a comparison of 
Figs. 33 and 34 (the ratio is closer to 2/3 here). For the 
purpose of this study, the kinematic omega fields are 
considered to best represent the actual vertical motion. 

A comparison of the vertical velocity fields with the 
horizontal mass transport shows consistent features. The 
vertical velocity maxima (Fig. 34a) correlate well tempo- 
rally with the surface convergence and upper level 
divergence maxima (Fig. 29a), and with periods of rapid 
surface deepening (Fig. 2b). In addition, the more intense 
vertical motion at radius 6 (Figs. 32, 33 and 34) , which is 
roughly 1.5 times greater than at radius 10, agrees well 
with the satellite discussion which depicted extreme amounts 
of deep convection in the inner radius (see Chapter 4) . 

E. CCNCLUSICNS 

In summary, stability and QLD mass budget analyses have 
revealed the following insights: 

• Surfaoe stability decreases during the early time peri- 
ods favor organized convection; 

• The stability analysis shows destabilization trends 
consistent with numerical model simulations (Sandgaths, 
1981) ; 
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• Evolution of cold core characzer isxics is svidenz 
throughout irost of the troposphere afzer the surface 
circulation intensifies; 



• Maxima in surface layer co 
ods of intense vertical mo 
divergence, coincident wi 
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sive development; 
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• Vertical velocity fields employing kinematic calcula- 
tion methods represent well rhe LND and vertical mass 
transport which is a vital link in the lifeblood of 
storm dynamics. 



• Vertical-radial cross sections of mass transport during 
extratrcpical cyclogenesis suggest similar development 
characteristics to tropical storm rapid intensirica- 
tion . 
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VI. 



VOfiT ICITY BU DGET ANAL YSIS 

A. GENERAL 

The absolute vorticity budget analyses are presented by 
vertical-time sections and by selected vertical profiles. 
Time periods refer to the 12-h span between synoptic times. 
For example, period 1206 refers to 0000-1200 GMT 12 January. 
In figures not referencing time periods, times refer to the 
specific time (GMT) . The 775 mb layer refers to the layer 
between the 850 and 700 mb. The sign convention is positive 
for processes producing vorticity increases in the budget 
volume (sources), and negative for vorticity decreases 
(sinks) . 

The main feature of this section includes a discussion 
of the absolute vorticity time tendency as a measure of 
storm development. Specifically, an investigation is con- 
ducted intc the role of horizontal and vertical absolute 
vorticity transport, advection, generation and dissipation 
in storm dynamics. Horizontal transport terms are separated 
into mean and eddy mode oomponents to gain better insight 
into their speoific contributions toward storm development. 
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Vertical transport terms reflect vertical divergence and 
advection frcm cmega estimates using the kinematic method. 
Referral to Section B of Table V is recommended as specific 
budget terms are identified and discussed. 

B. TIME SECTICNS , TIME TENDENCIES AND SELECTED PROFILES 

The dramatic buildup cf absolute vorticity in the budget 
volume at both radii over the 2-day span of intense 
development is illustrated in Fig. 35. At 0000 GMT 13 Janu- 
ary (1300 GMT), the relative cyclonic vorticity in the inner 
volume through the 600 mb level is not much different from 
the ccriolis farametar (7.82 X 10**-5). Notice the presence 
of anticyclonic relative vorticity aloft, which is associ- 
ated with the weak and dissipating upper level ridge. The 
dramatic increase in absolute vorticity by 0000 GMT 15 Janu- 
ary in the inner radius occurs from two sources. The rapid 
movement to the northeast over the period has resulted in a 
larger value cf planetary vorticity (f = 11.56 X10**-5), 

while the increased storm circulation accounts for the 
additional buildup. Although the surface circulation is 
strong, the naximum absolute vorticity exists in the middle 
troposphere. Similar changes occur at radius 10, although 
the magnitude cf the absolute vorticity increase is smaller. 
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and there is a uniform increase through the ^.roposphere 
below the jet core level. Specific factors contributing to 
the worticity changes at the different vertical levels will 
be discussed as the specific budget terms are analyzed. 

A different perspective of the absolute vorticity 
buildup in the budget volume is presented in Fig. 36. Here, 
vertical time sections of vorticity for both radii illus- 
trate several key features. At radius 6 (Fig. 36a), a vor- 
ticity minimuiE can be seen at 200 mb during the early 
period. The incr^se with time at lower levels is a reflec- 
tion of the nascent cyclone vertical structure which is com- 
posed of anticyclonic vorticity aloft over the low-level 
developing system. After 1212 GMT, absolute vorticity 
increases throughout the whole troposphere until 1500 GMT. 
Two significant features are evident during these rimes. A 
large buildup occurs in the upper troposphere (vicinity 225 
mb) after 1300 GMT, which causes a temporary reduction in 
the vertical vorticity gradient. A distinct lower tropos- 
pheric maximum is evident at 1500 GMT, which occurs just 
after the period of intense circulation increase, and is 
substantiated by the mass transport maxima shown in the mass 
budget analysis (Fig. 29a) and the associated central 
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surface pressure falls illustrated in Fig. 2a of the synop- 
tic discussion. This large maximum in vorticity is coinci- 
dent with the period of lowest central pressure and maximum 
surface layer circulation. Vorticity values are seen to 
decrease in tfce middle and lower troposphere for 24 h after 
1500 GMT, but at a significantly slower rate than the 
buildup. Values increase again through the whole 

troposphere after 1600 GMT due to the inclusion of the 
developing system to the south. 

Several important differences exist at radius 10 (Fig. 
36b). Vorticity values increase at a slower rate and at 
later times than in the inner volume. This result indicates 
that although vorticity values increase in both volumes dur- 
ing the development stage, the outer volume experiences a 
less significant increase due to the greater spatial extent 
and relatively smaller vorticity. 

Absolute vorticity time tendency, which is the first 
term in the vorticity budget equation (see Table V) , is pre- 
sented in Fig. 37. This display is a time derivative of the 
fields depicted in Fig. 36. Recall that time periods, when 
used in subsequent figures, refer to the 12-h time interval 
between synoptic times. 
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A large increase in the vorticity tendency in the inner 
volume occurs in the upper levels from 1218-1318 GMT. Vor- 
ticity increases at a constant rate in the surface layer 
during the development stage. The general decrease in ten- 
dency values through the whole troposphere after 1506 GMT is 
consistent with the onset of the occlusion phase and storm 
dissipation. It is interesting that a vorticity tendency 
minimum occurs near 450 mh during a period of large surface 
pressure falls (Fig. 2b) and maximum mass transport (Fig. 
29). No correlation to synoptic causes for this tendency 
minimum is apparent from detailed inspection of the 500 mb 
height and vorticity fields (Figs. 4b, 7b, 10b). One inter- 
pretation is that a subsynoptic scale system in the wind 
fields affects the vorticity values at these levels and 
times. An alternative cause could be linked to insufficient 
data input into the 1312 or 1400 GMT analyses. 

Similar significant features exist at radius 10 (Fig. 
37b) . The nature of the buildup of absolute vorticity 
suggests inherently similar mechanisms are responsible for 
the absolute vcrticity buildup between inner and outer vol- 
umes. Mechanisms such as the jet core influence, mid-trc- 
pospheric self-amplif icaticn, and horizontal and vertical 
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•transport are possibilities explored in the following 
sections. 

C. LATERAL TRANSPORT 

The lateral transport of absolute vorticity is one of 
the three most significant forcing terms in the budget equa- 
tion (Table V) . Transport into the budget volume arises 
from the mean convergent flow (mean mode) and from 
asymmetries in the cyclonic flow (eddy mode) . 

A typical view of the lateral vorticity transport occur- 
ring during storm development is presented in Fig. 38. The 
curves are valid during the maximum development period at 
radius six. The mean mode and total lateral transport 
curves are similar and show similar trends. The mean mode 
is seen to dominate the total transport at levels below 800 
mb. In the upper levels, the total lateral transport is 
small, but is the residual between two strong oppositely 
contributing processes. These features constitu-ce signifi- 
cant differences from those observed by Conant (1982) in the 
President's Day storm, where he observed than the total and 
eddy mode lateral transport curves were the more strongly 
correlated. Although the curves presented in that study 
represented an earlier phase in the storm development, the 
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general features indicate that different mechanisms serve -o 
increase inner volume vorticity levels within each storm 
budget volume. A look at the mean and eddy mods lateral 
transport vertical-time sections gives insights into 
processes responsible for these differences. 

Another interesting feature is the presence of a measur- 
able eddy mcde contribution in the lower troposphere, which 
is an indication that processes other than mean convergence 
are acting to increase the vorticity and storm circulation 
within the six degree volume. These results should be 
expected, since processes such as intense convection in pre- 
ferred quadrants of the storm, frontal asymmenr'ies and 
elliptical vortex structure may cause asymmetric flow about 
and into the cyclone. Thus, some eddy mode contribution is 
evident through the whole troposphere. Mean and eddy mode 
contributions to the total vorticity lateral transport are 
discussed below. 

The total lateral vorticity transport time sections for 
the inner and outer budget volumes are presented in Fig. 
39. In the inner volume (Fig. 39a), inward transport exists 
below 775 mb through all time periods, with alternating 
inward and outward transport evident in the upper 
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troposphere (near 350 mb) . Notice the mid-tropospheric (600 
mb) outward transport (after 1506 GMT) as broad divergence 
is established. The outward transport maximum at 225 mb 
suggests strcng divergence aloft as the jet stream dominates 
the upper troposphere flow. The surface layer regime is 
particularly noteworthy in that especially deep inward vor- 
ticity transport occurs during the deepening phase 
(1318-1506 GMT) through U50 mb and contains a maximum (19.2 
X 10**- 10/sec**2) at 925 mb. This signature is entirely 
different from that observed by Conant (1982) in the Presi- 
dent's Day storm. He showed a maximum at 220 mb with rela- 
tively weak surface transport. The emphasis on the level of 
maximum horizontal lateral transport serves to focus on the 
differences between the underlying processes responsible fcr 
vorticity increases within the budgets of both systems. In 
this case, the lateral transport surface maximum occurs in 
temporal coincidence with the largest central pressure falls 
and maximum mass transport, as discussed in previous 
chapters. 

The total lateral transport at radius 10 (Fig. 39b) 
exhibits similar inward transport in the surface layer 
although the magnitude is considerably less and the maximum 
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occurs a time period earlier. Inward vorricity transport is 
evident through the rest cf the troposphere (except abcve 
450 mb after 1606 GMT ) , with a secondary maximum observed 
between 450 and 175 mb during the 1218-1318 GMT periods. 
The time and level of occurrence of this maximum strongly 
suggests vcrticity import from jet streak interaction. Sub- 
sequent low-level inward transport occurs from circulation 
increases in response tc this upper-level forcing. The 
earlier occurrence of both upper and lower-level maxima in 
the cuter volume suggests earlier forcing from jet streak 
interaction as it enters the outer volume first and subse- 
quently migrates toward the inner volume. The increase in 
the surface layer transport at 1618 GMT is again due to the 
developing system to the scuth. 

D. MEAN MCDE LATERAL TRANSPORT 

As described in Chapter 3 (Table V) , the absolute vor- 
ticity transport can be partitioned into mean and eddy 
modes. The purpose of this partitioning serves tc isolate 
contributions by the mean flow (mean cyclone convergence and 
divergence), and the perturbation flow (the covariance of 
vorticity and normal wind component deviations) . 
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I. 





Specifically, the asan cyclonic low-level inflow and upper- 
level outflow of vorticity is represented by the mean mode 
and reflects necessarily the convergant/divergent nature of 
the Icwer/upper layers. 

At radius six (Fig. 40a) , the strong dual-centered 
transport maximum, the upper-level outflow maxima and the 
level of transition between inward and outward transport 
agree well both temporally and spatially with the mass 
budget transport features discussed in Chapter 5. Periods 
of maximum surface cyclone deepening coincide directly with 
the vorticity transport maxima discussed above. Differences 
in mean mode transport between the President's Day storm and 
the present case arise due to the significantly weaker 
surface-layer lateral vorticity transport in the President's 
Day storm, whereas the upper layer transports are nearly the 
same. Thus, the main difference in the growth and vitality 
in this storm is the significantly higher mean flow 
transport of vorticity in the lower layers. 

At radius 10 (Fig. 40b), absolute vorticity lateral 
transport values are relatively weaker than at radius 6. A 
single inward transport maximum in the surface layer and a 
coincident upper-level outward maximum are evident. These 
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in response to earlier dynamic 



maxima occur earlier, again, 
forcing in the outer volume by the upper-level eddy 
structure changes (subsequently addressed) . The effects of 
the new lew to the south are again seen after 1518 GMT. 

E. EDDY MODE LATERAL TRANSPORT 

The eddy mode absolute vorticity lateral transport gives 
insight into the upper level vorticity increases which ulti- 
mately affect the lew level circulation tendencies. As 
discussed above, the eddy mode component represents the hor- 
izontal transport due to asymmetries in the cyclonic flow 
induced by wind field correlated vorticity deviations. 

Inspection of the eddy mode absolute vorticity transport 
for radius 6 (Eig. 41a) shows considerable vorticity import 
in the layer between 750 and 125 mb after 1206 GMT. Within 
this broad area of positive lateral vorticity transport are 
two significant time periods. Between 1206-1306 GMT, a max- 
imum of vorticity transport is seen to exist in the 450-275 
mb layer. Inspection of the 300 and 400 mb vorticity and 
winds for 1200 and 1212 GMT shows weak ?VA into the budget 
volume, and this is substantiated by the coincident maximum 
in the calculated advection (Fig- 43a) discussed subse- 
quently. These facts suggest vorticity import in the form 
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of the short wave that tecomas established over southern 



Japan during this time frame (discussed earlier in the 
synoptic discussion. Fig. 4b). The effect of this perturba- 
tion diminishes with time as the cyclone accelerates to the 
northeast: away from the slower moving upper-level trough 
after 1300 GMT, and thus, the eddy mode transport also 
diminishes. Vorticity increases are again evident after 
1306 GMT when the jet streak has become well established in 
the budget volume (Fig. 7a). The influence of the jet 
streak continues until 1518 GMT when the jet core becomes 
detached from the budget volume as the surface center 

migrates steadily northward. 

Significant eddy mode contributions are restricted to 
the upper levels, and are only roughly 1/2 those observed by 
Conant (1 982) in the President’s Day storin. This is an 
extremely valuable result, as Conant showed that the migra- 
tory short wave played a major role in that szcrm 

development. In this study, the migratory short wave is 
shown to be of minimal influence, and thus vorticity 

increases in the upper levels occur primarily from cyclonic 
shear vorticity input into the budget volume as the jet core 
becomes established. Shile the eddy mode contributes 
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considerably to the initial storm development, this effect 
is net nearly as dominant during the explosive stages as in 
the President's Eay storm. In this case then, it is logical 
to conclude that explosive development occurs due to upper- 
level jet forcing acting in concert with the surface 
development along the barcclinic zone. The secondary eddy 
mode maximum near 925 mb at 1418 GMT is presumed to be due 
to asymmetries in the mean flow in specific quadrants of the 
storm circulation. 

Similar features exist at radius 10 (Fig. 41b). As in 
the inner volume, significant contriburions occur az the jet 
core level with maxima of similar magnitude occurring at 
coincident times. A weak maximum exists from 1406-1506 GMT 
in the surface layer, and as in the 6-degrae volume, 
represents small asymmetries in the mean cyclonic flow. 

Alternatively, the lateral transport may be expressed in 
terms of lateral divergence and advection components through 
the use of vector identitites (see Table V) . These terms 
provide a mote traditional view directly in terms of the 
vorticity equation. The divergence component of the parti- 
tioned lateral vorticity transport (Fig. 42a), is composed 
of a two-layer regime with positive vorticity production due 
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to convergence in the surface layers and negative vorticity 
production with divergence aloft. These features match 
closely the associated mass transport features discussed in 
Chapter 5. Ihe temporal and spatial occurrence of the 
inflow and outflow meixima agree well with similar features 
of the mean mode discussion, and with the general 
circulation increases and maximum pressure falls. 

The most interesting comparison occurs between the eddy 
mode and horizontal advection transports (Figs. 4la and b 
and 43a and b) . The upper tropospheric values match well 
both temporally and spatially, as well as in magnitude. 
However, the advection term in the lower troposphere shows a 
substantial negative vorticity advection (NVA) into the vol- 
ume below 775 mb, which is quite different from the negative 
and positive eddy vorticity transports. In-depth examina- 
tion of the vorticity and wind fields in the inner volume 
during a period of maximum development shows an asymmetric 
east-wesr oriented vorticity field coupled to a strong 
cyclonic inflow pattern, which causes a net advection of 
smaller vorticity values into the volume in the northeast 
and southwest quadrants. This net NVA is further enhanced 
by the 40-kt translation speed of the surface center to the 
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northeast during this period. The change zo a westerly, 
more zonal wind component with height accounts for the sign 
reversal in the vorticity values above 775 mb. Notice that 
the addition of the advection and divergence terms is very 
close to the sum of the mean and eddy mode lareral 
transports. These results serve to verify the integrity of 
the transport calculations. 

In summary, the total horizontal absolute vorticity 
transport has teen shown tc be composed of a mean mode rep- 
resenting contributions from the mean mass circulation, and 
an eddy mode representing contributions from asymmetries in 
the flow from advection and jet core sources. The mean mode 
transport was shown to be very similar to the total 
transport through all levels, while the eddy mods transport 
was shown to be significant only above 500 mb. The rela- 
tively smaller magnitude of the eddy mode transport: maxima 
compared to thar observed in the President's Day storm 
suggests this term plays a less significant role in the 
storm evolution. 

F. VERTICAL EEDISTRIEOTICN 

This section investigates the role of vertical vorticity 
transport from the upper and lower layers into the middle 
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troposphere (as discussed earlier in connection 
34) . In this case the diffluent nature of r 
evolving trcuch will be explained in terms of this 
vorticity increase. The concentration of vortic 
mid-troposphere serves to increase mid-level ci 
increase theriral advection, and amplify th 
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The divergent component (Fig. 45a) is the main 
contributor tc vertical vcrticity redistribution. This term 
is simply the opposite of the lateral divergence term (Fig. 
42a) , but differs slightly due ro the kinemaxic vertical 
motion estimates. There is a decrease of vorticity below 
the maximum in the vertical motion, and a gain above this 
maximum. The reversal in sign from the lateral divergence 
term (Fig. 42) is consistent from continuity, and close 
agreement between the tempcral and spatial magnitudes shows 
consistency. Similar features exist at radius 10, although 
the magnitudes are considerably less as expected. 

The vertical advection portion (Fig. 46a) plays a rela- 
tively minor role in the budget. Positive vertical advec- 
tion is seen tc occur between 775 mb and 350 mb from 
1318-1506 GMT. Accumulated values of surface layer vortic- 
ity are advected by the vertical motion during these times. 
Positive values below and negative values above 450 mb con- 
tribute to the vorticity buildup in the middle troposphere 
illustrated in Fig. 35. The relatively minor contribution 
of vertical advection to the budget results is consistent 
with results found by Conant (1982). 
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In summary. 



the vertical 



transport redistributes 



large vorticity buildup in the 
mid and upper troposphere. The 
major contributor, as the advect 
throughout the budget volume, 
similar with considerably lower 
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An impcrtant feature tc notice is the relatively large 
magnitude of the maxima in both regimes compared tc the 
total lateral transport. These patterns suggest -hat the 
divergence term plays the dominant role in the budget. In 
the lower troposphere, convergence overwhelms all other 
processes, as is expected from a scale analysis of the 
general vorticity equation. 

The tilting term (Fig. 47a and b) arises from vertical 
vorticity components generated by the tilting of horizon- 
tally oriented vorticity elements by a non-uniform vertical 
motion field. Strong vertical shear and a horizontally var- 
ying omega field are prerequisites for contributions from 
this source. It is apparent (from Fig. 47a and b) that the 
products of vertical shear and the change in the vertical 
motion fields in the horizontal are small compared to 
lateral and vertical transports. 

Frictional dissipation (Fig. 48) is assumed to occur 
only in the surface boundary layer and is parameterized 
using a stability dependent scheme (Johnson and Downey, 
1976). The calculated values were about 1.5 times greater 
than those observed by Ccnant (1982), as is expected from 
the larger intensity cf the disturbance studied here. The 
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maximum at 1506 GMT has a magnitude very close to the value 
of the low-level maximum in the total lateral transport 
(Fig. 39a). These results suggest that frictional effects 
assume a prominant role in the dissipation of vorticity in 
the surface layer. Values in the outer volume (Fig. 48b) 
suggest minimal frictional influences on the vorticity 
averaged over a larger region. 

In summary, the divergence term was a major factor in 
the cverall budget as a source of cyclonic vorticity in the 
lower troposphere and a sink aloft. The friction term was 
an important sink of vorticity in the surface layer, whereas 
the tilting term had no significant effect on the budget. 

H. EESIDOALS 

The residual in the vorticity budget contains the 
effects of the coitted processes plus the accumulated errors 
in the calculation of the resolved terms. Inaccuracies in 
the vertical mcticn and horizontal wind fields contribute to 
the physical residual components, whereas spatial and tempo- 
ral finite differencing and grid-point interpolation inaccu- 
racies are causes for ccmputat ional errors. A positive 
vorticity residual indicates an apparent vorticity source in 
which the observed vorticity increases at a point are larger 
than is estimated by the computed terms. 
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I. CONCLUSICNS 



In 


summary 


, the 


increase of absolute 


vorticity in 


the 


budget 


volume 


during 


th <= 


development 


St 


age was shown 


to 


occur 


from a 


complicated 


interaction 


of 


process es . 


The 



period of surface layer circulation organization associated 
with incipient development is distinctly different from the 
period of explosive development in which the circulation 
increases throughout the troposphere. 

During the incipient stage (period 1206-1318), the lower 
troposphere circulation increases were shown to occur prima- 
rily from the source divergence term, as the storm organized 
along the fcarcclinic zone. Lateral transport increased as 
the circulation developed and contributed to small net vcr- 
ticity increases since the frictional dissipation was also 
enhanced. 

Mid-troposphere vcrticity values increased in the earli- 
est period due to weak lateral transport in the form of PVA 
from the short wave established over Japan (discussed in 
detail in earlier sections) . This mid-level vorticity 
import acted to enhance the low-level organization in the 
earliest stages. At later times prior to 1318 GMT, the 
divergence of the vertical transport was seen to account for 
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a significant portion of the mid-level vorticity increases. 
An increase in the deep divergence layer above the LND 
served to offset a portion of the mid-level increases. 

In the upper troposphere, vorticity values increased as 
the jet core propagated into the budget volume. Additional 
increases of snaller magnitudes were evident from vertical 
redistribution from the surface layer. These two processes 
combined to more than offset negative vorticity due ~o the 
increased divergence. Thus, upper level forcing was 
available for surface development. 

The most significant changes to account for the explo- 
sive development after 1318 GMT occurred in the upper 
troposphere. The small positive lateral transport was a 
balance between the mean mode (divergence) and eddy mode 
(advection) contributions. Large amounts of cyclonic shear 
vorticity were advected into the budget volume after 1306 
GMT as the jet streak migrated into the budget volume. The 
corresponding divergence increases occurred 12 h later, 
which served to promote the explosive low-level development. 
The maximum divergence established a time period later 
served to ventilate the system and induce the maximum mass 
outflow aloft, increasing the surface circulation, and 
promoting the explosive development. 
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Kid-tropo sphere 



vorticit y 



increases 



after 



1318 GMT 



resulted priiarily from the divergence of rhe vertical 
transport and the upward development of the mid- level 
diffluent trough from thermal advection. The satellite 
sequence in Chapter U indicates that convection becomes 
intense and widespread during the period due to vertical 
stretching in the mid-trcposphere induced by the increased 
upper-level divergence. Vertical stretching acts zo cool 
the middle troposphere and destabilize the vertical columns. 

In the lower troposphere, intense circulation increases 
resulted from the induced upper-level forcing associated 
with the jet streak trovement into the budget volume. The 
net vorticity increases results from the excess of the 
lateral transport over the frictional dissipation. Although 
the frictional effects increased at about the same rate and 
to the same magnitude as the lateral transport, there was a 
lag cf 12 h. When the upper level forcing and ventilation 
decreased as the surface center was removed from the jet 
influence, the system dissipated after 1506 GMT. 
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VII 



CONCLUSICNS and RECCHMEND ATION S 



A synoptic investigation was conducrad of a Western 
Pacific case of extratropical explosive cyclogenssis during 
12-17 January 1979. Storm related mass and vcrticity budg- 
ets were computed for the cyclone using the SCMWF FGGE 
analyses. Overall storm development features were examined 
and discussed. 

Ihe following significant results ware developed; 



• Maximum circulation increases due to rapid pressure 
falls match well temporally and spatially to maxima in 
the surface layer inward mass transport (convergence) 
and upper-level outward transport (divergence). 



• PVA in the mid- troposphere in the form of short wave 
interaction, in combination with thermal advection 
effects, enhances the organization of The low-level 
circulation in the initial stages of development; 



• The initial short wave interaction decreases after the 
initial period while surface development continues. 
Explosive development occurs coincidentally with the 
migration of a strong jet streak into the budget vol- 
ume, which Drovides the necessary upper-level forcing 
from divergence and shear vorticity advection. The 
lack of an upper-level migratory short wave evident 
during explosive development confirms Petterssen type-A 
development ; 



• The stability analysis reveals trends consistent with 
numerical simulations (Sandgathe, 1981). Stability 
decreases rapidly during surface development as cold 
air was drawn over the warm sea surface. Stability 
values increased after maximum storm development, as 
cold air advection dominated the lower troposphere. 



• Mass transport radial cross sections suggest explosive 
development characteristics similar to those observed 
in tropical storm intensification with the followina 
similarities: 1. low static stability in the surface 

layer, 2. the availability of an oceanic heat source, 
3. strong outflow (divergence) aloft, and 4. strong 
convection rear the inner radii; 
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• The FGGE data set was found to be Isss than ootimum for 
supporting 6-h QLD budaet analyses. The 060D and 1800 
GMT time periods contain spurious potential temperature 
biases affecting all of the budget calculations. The 
problem alsc affects both initialized and uninitalized 
fields which sugaests inaccurate analysis inputs to the 
model first guess. 



Eecommendations for further research include: 



• The O'Brien correction scheme applied to the horizontal 
mass transport should be examined to establish a more 
optimum representation of the vertical motion fields; 



• Further efforts should be made to fully assess the 
impact of the FGGE 0600 and 1800 GUT data problems and 
their effects on represented fields in the FGGE data 
set. 



• Energy diagnostic budget calculations should be ccm- 

f leted on fhis cyclcne to isolate the contribution of 
hermal advection, sensible heat and latent heat in the 
explosive development. Correct parameterization of 
these inputs could serve as a basis for further compar- 
isons between explosive extratropical and tropical 
storm formation. 
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APPENDIX A 



FGGE DATA EIAS CHARACTERISTICS 

The fcllcwing information is presented to document the 
6-h FGGE data problem encountered in this study. As indi- 
cated in Chapter 5, a significant high potential temperature 
bias was cbserved in the 850-500 mb layer at 1800 GMT 14 
January. Several additional fluctuations were observed to 
occur at 0600 and 1800 GMT on other days. The nature of 
these temperature biases radically affected the budget cal- 
culations and precluded the planned use of a 6-h interval 
format. The nature of the effects is shown in Figs. 50 and 
51. The top portion of Fig. 50 is a 6-h representation of 
the vorticity time tendency term at radius six, whereas the 
bottom field is the inner volume budget residual. The sig- 
nificant 6-h oscillations in both fields are indicative of 
the difficulry encountered using QLD budget techniques on 
this particular FGGE data set. 

An even more striking example of the complexity of the 
6-h data problem is shown in the inner volume stability 
trace (Fig. 51). Notice the oscillation between maxima and 
minima indicating the alternate warming and cooling 
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tendencies in the 6-h fields. It is obvious that a synoptic 
time scale interpretation of this figure is virtually impos- 
sible. The 1300-1500 GMT period is especially interesting. 
The 0600 and 1800 GMT stability values are markedly larger 
than at 0000 and 1200 GMT. The increased stabiltiy during 
these times occurs from warm biases in the mid-level 
potential temperatures, as shown in Fig. 27a and b. 

To determine possible causes for this 6-h data problem, 
specific mid-tropospheric observational data available as 
input to the FGGE analysis scheme at 1418 GMT was studied. 
The observations consisted solely of satellite input (micro- 
wave and infrared wavelengths). Close scrutiny of the sat- 
ellite input in ccmpaiison with the FGGE gridpoint estimates 
revealed interesting discrepancies. The final FGGE 

potential temperature analysis was 5-7 K higher than the 
associated satellite data in the southwest quadrant of the 
budget volume. In all other regions the soundings and the 
FGGE fields matched vary closely. In addition, potential 
temperature biases of smaller magnitudes were evident at the 
1800 GMT time period on ether days (Fig. 27) . 

Further insight is gained by an inspection of the unini- 
tialized height fields (through the thickness analyses in 
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Figs. 16c, 19c and 22c) . These fields also show flucrua- 

fions of the same scale occurring at the 1800 GMT time 
frame, as this thickness field departs markedly from the 
preceding and subsequent analyses. It is obvious rhat the 
6-h data problem affects both initialized (potential temper- 
ature) and uninitialized (height) fields. The evidence pre- 
sented above suggests that a problem existed in the analyzed 
potential temperature fields used in the first guess, and 
that the nature cf the problem was not identified during the 
FGGE data post-processing evaluations. 

Further support of these conclusions is presented by 
Paegle (1983), who summarized features of the FGGE data set 
in tropical regions. Lower tropospheric temperature biases 
on the order cf 2-3 degrees were found for some time peri- 
ods, which were presumably from the lack of conventional 
data input at the 1800 GMT times. In addition, the SCMWF 
forecast model was shown to be insensitive to diurnal temp- 
erature changes over land. To what degree this is evident 
over oceanic areas is still uncertain. 

In summary, it is believed that the spurious warming at 
1800 GMT between 850 and 500 mb is due to a problem with the 
FGGE analysis scheme rather than a lack of quality 
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obssrvatiocs. Ihe inpact of this assessment on the resul-s 
of this thesis cannot be fully ascertained from the prelimi- 
nary cbservaticns discussed here. Detailed investigation 
into specific analysis fields used in the first guess must 
be dene level by level tc fully assess the impact of this 
6-h data problem. 
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APPENDIX B 



TABLES 
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20 mb 
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(From Eengsston, 1982) 
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TABLE II 

Lis-ir.g of Analyses and Prognoses Utilized 



(FGGE=1 ,NMC=2, PN0C=3) 



LEVEL 


HU 


ANAL/PROG 


ORIGIN 


-level 


pressure 


anal 


1,2,3 


-level 


pressure 


36-h prog 


2,3 


-surface 


te Diperature 


anal 


3 


-mb 


he ights 


anal 


1,2,3 
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TABLE III 

Generalized QLD Budget Equation in Isobaric Coordinates 

(After Wash, 1978) 



1000 nb 0 0 

where F is the volune integral of the desired budget property f. 
The budget equation is 

= LT(f) + VT(f) + S(f) 

where the lateral transport is 
100 nb 2:: 



100 nb 6 2 tt 




LT(f) 



(U-W)g frsinB da(-dp)l^ 



1000 nb 0 



and the vertical redistribution is 



100 nb 3 2 it 



VT(f) 



/ 




1000 nb 0 0 



The source/sink tern is 



100 nb 3 2ir 




1000 nb 0 0 
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TABLE IV 



QLD Mass Budget Equations in Isobaric Coordinates 

(After Wash, 1978) 

The Definition 



■■/i 



r sinB dadSdp 



where f = 1 . 



The Budget Equation 



^ = LT + VT , 
at 



T tV> /-»•>» /n 



Lateral Transport 



h™ 2 tt 



LT = - 




(U-W) ^ r sinS^ dadp ] g 






Vertical Transport 



3^ 2r 



VT = 




1 2 

— (oj-(o„) r sinB dadB ! 
g B P 



0 0 



where 



dp 
dt ’ 



dp. 



= 



B dt 



96 



TABLE V 



QLC Circulation Budget Equations in Isobaric Coordinates 

(After Wash, 1978) 



Section A 



C = 
a 



100 mb B 2 tt 

/ II i 



C r sin3 dadS (-dp) 

3. 



1000 mb 0 0 



where C is the absolute circulation and ? is absolute vorticity. 

3 - 3 



The budget equation is 



£C 



= LT(? ) + DVT(;^) + S(; J 

O L 3 3 3 



where the lateral transport is 
100 mb 8 2 tt 

LT(?a) = y* / / 6 da(-dp)lg 

1000 mb 0 0 

and the divergence of the vertical transport is 
B 2tt 



■// t 



D\n:(^ ) = / f (w c ) r^ sinB dadB . 

a II g dp a 



0 0 



The source/sink term is 
100 mb B 21T 



S(?^) 

a 



■I II i 

1000 mb 0 0 



dC 



— sinB dadB (-dp) 
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TABLE V (Cont.) 



Section B 



The partitioned form of the vorticity budget equation is 



6(C) 

a 

at 



+ LT (5 ) + DVT (? ) 

a , a 



+ s 

a 



1 

mode 


1 

f vertical 

eddy mode divergence 


r 

vertical 

advection 




f" ' 


! 


1 


1 



horizontal horizontal 
divergence advection 



divergence tilting frictional 
term term dissipation 



The above partitions make use of Stokes’ theorem 



/ 



? U • m dZ = 
a ~ ~ 






dA , 



and the division of total flux (U ^ ) into divergent and advective 
components , 



U = S (V*U) + U*V S 
a ~ a ~ -- a 



98 



TABLE VI 



Mean and Eddy Bode Transporx Equations 
(After Conant, 1982) 






a 



EM(g^) 

3 . 



3 



100 nb 2 tt 

1 ■"* 

~ ^ da(-dp)|g 

1000 mb 0 
100 mb 2^^ 

~ ^ (U-W)“ r sinB da(-dp)L 
g a - ~ p P 

1000 mb 0 
100 mb 2 t 

1 “* 

~ ^ da(-dp)lg 

1000 mb 0 
100 mb 2ir 

~ (U-W)p r sinS da(-dp)L 
g a - - P p 

1000 mb 0 







EM is the eddy mode and MM is the mean mode. 

2tt 



GC 

Note that ( ) 




( ) da 



and ( ) = ( 

from its mean around 



/ — ■ ct 

) - ( > is the deviation of the property 

a lateral boundary. 
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SFPENDIX C 



FIGURES 




gure 1. 



Th€ Overall Data Hanagement During FGGE (From 
Bsngsston, 1982). 
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Figure 2. 



Stcrm Track Positions and SLP Analyses. (A) 
Storm Positions in 12-h Increments. (B) Sea 
Level Pressure Analyses in 12-h Increments. 
12/0000 and 1200 refer to 0000 GUT 12 Januar 
1979. 



Times 

y 
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Fiqure 3. 36-h Forecast SLP and Storm Position Errors. (A) 

S-a Level Pressure Forecasts in 12-h Increments- 
(E) Position Errors Based on FGGE Positions. 

Time 1200 refers to 0000 GMT 12 January 1979. 
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Figure 4. Svncptic Fields for 0000 GMT 12 January 1979. 

(A) winds/Isot achs (m/sec) . (B) 500 mb Absolute 

vcrticity (DASHED) X 10**5/sac and Heights 
(SOLID) rn gpm. (C) Sea Level Pressure (SOLID) in 
mb and 1000-500 mb Thickness (DASHED) in gpm. 
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Figure 5. 



DHSE Visual Satellite Imaaery for 0041 GMT 12 
January 1979. labels are ?A) Incipient System 






Secondary Cevelcping 
“ream. 



ys' 

ystem (C) PVA Area 
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Figure 6. 



DaSE Infrared Satellite Imagery for 0041 GMT 12 
January 1979. labels are (A) Incipient Svsrem 
(E) Secondary Developing 
upstream. 



ystem (^) PVA Area 
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Figure 7. Same as Figure 4 except for 0000 GMT 13 January 
1979. 
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Figure 8. 



Same as Figure 5 except fo 
1979. Label A is Incipient 



2341 GMT 12 January 
System . 
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Figure 9. 



Sams as Figure 6 
1979. Label A is 



except for 2341 GMT 
Incipient System. 



12 January 
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Fiqura 10. Same as Figure 4 except for 0000 GMT 14 January 
1979. 
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FNCC SST Analysis Valid 1200 GMT 14 Janaary 
1S79. "X" marks cyclone position at 0000 GMT 

January. Isotherms in degrees Celsius. 
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Figure 11. 




Figure 12. Same as Figure 5 except for 1405 GMT 13 January 
1S79. Label A Indicates Jer-Frontal 
Intersection. 
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Figure 13. Same as Figure 6 except for 1405 GMT 13 January 
1979. Label A Indicates Jet-Frontal 
Intersection . 
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Figure 14. 



Same as Figure 5 except for 0023 GMT 14 January 
1S79. Label A Indicates Feeder Bands. 
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Figure 15. 



Same 

1979 



as Figure 
Label A 



6 except for 0023 GMT 14 
Indicates Feeder Bands. 



January 
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Figure 16. Same as Figure 4 except for 1200 GMT 14 January 
1 S ^ S • 
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Same as Figure 5 except for 1205 GKT 14 January 



1 16 



Figure 17. 







Figure 18. 



Same as Figure 6 except for 
1S79. 



1205 GMT 14 January 
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Fiqure 19. Same as Figure 4 except for 1800 GMT 14 January 
^ 1979. 
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Figure 20. Figure 5 except 



1 19 



for 2124 GMT 14 January 





Figure 21. 



Same as 
1S79. 



Figure 6 except for 2124 GMT 



14 January 
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Figure 22. 



Same as Figure 4 except 
1S79. 



for 0000 



GMT 15 January 
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Figure 23. 



'^ame as Figure 4 axcepr for 0000 GMT 
1S79. 



16 January 
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Figure 24. Same as Figure 5 except for 1005 GMT 16 Januarv 

I V / M 
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Figure 25. 



Same as Figure 6 except 
1S79. 



for IOOd GKT 16 January 
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Fiqure 26. Same as Figure U except for 0000 GMT 17 January 
^ 1979. 
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Figure 27. Area-Averaged Potential Temperature Time 

Sections (A) Radius 6 (B) Radius 10. Values are 
in Degrees Kelvin. Time 1200 Refers to 0000 GMT 
12 January. 
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TIME(GMT) 




Figure 28. Stability Index (A) Observed (B) Simulated (From 
Sandgathe, 1981). Values in Degrees Kelvin / 

100 mb. Time 1200 Refers to 0000 GMT 12 
January. 
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Figure 29. Ccrrected Lateral Mass Transport Time Sections 

(a) Radius 6 (E) Radius 10. Contour Interval is 

40 X (10 * ♦ 'll). Dashed (Negative) Values 
Indicate Flux Cut of the Volume. Units are in 
am/sec-100 mb. Period 1206 Refers to 0000-1200 
Out 12 January. 
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30, 



Radial Section of Lateral 
Time 1400 ]B) Time 1412* 
in Figure 29. 



Mass Transport (A) 
Units and Contours 



as 



129 



125 

175 

225 

275 

350 

^0 

600 

775 

925 

src 

(25 

(75 

225 

275 

350 

450 

600 

775 

925 

src 

6 . 





1, Radial Sectxon 
Time 1500 (B) 
in Figure 29. 



of Lateral 
Time 1512. 



Mass Transport (A) 
Units and Contours as 
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1218 1306 1318 H06 1418 1506 






1518 1606 1618 




32* QLD Derived Vertical Velocities (A) Radius 6 (3) 

Radius 10, Units are in mb/lOOO-sec, Time 
Periods as in Figure 29, 
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Figure 33. 



FGGE Derived Vertical Velocities (A) Radius 6 
(E) Radius 10. Units are in mb/1 OOO-sec . Time 
Periods as in Figure 29. 
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Figure 34. Kirematic Vertical Velocities fA) Radius 6 (B) 

Radius 10. Units are in mb/1000-sec. Time 
Periods as in Figure 29. 
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Figure 35. Budget Volume Absolute Vorricizy for Hadiu 
and Radius 10 for Two Times. Time 1300 Z R 
to 0000 GMT 13 January. 



1 34 



W O' 




Figure 36. Area-averaged Absolute Vcrticity Vertical Time 
Sections for (A) Radius 6 and (B) Radius 10. 
Time 1200 Refers to 0000 GMT 12 January. Units 
are 10**-5/sec. 
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Absolute Vorticitv Vertical Tine-Tendency 
Sect'ons for (A) Radius 6 and (B) Radius 10. 
iericd 1206 Refers to 0000-1200 GMT 12 January. 
Units are 10 * ♦ -10/.(sec *.* 2). SoliQ 
Contours reflect vorticity increases. 
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Figure 38. Total, Mean and Eddy Mode Lateral Vortici 
Transport During Maximuni Develoomenr . T 
in Figure 35. unixs as in Figure 36. 
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Same as Figure 37 except for Lateral Vorticiry 
Transport . 
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Figure 40. 



Same as Figure 37 except for 
Transport. 



Mean Mode Lateral 
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!• Same as Figure 37 except for Eddy Mode Lateral 
Transport. 
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Figure 42. 



Same as Figure 37 except 
Divergence. 



for 



Lateral 



V ortici 



y 
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Figure 43 



Same as Figure 37 
Ad vection. 



for Lateral Vorticity 





Figure 44. Divergence of the Verrical Transport. Units, 
times and contours as rn Figure 37. 
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Same as Figure 37 except for Vertical Vorticity 
Divergence. 



144 




I2S-1 



(R) 



I7S- 



02S- 



SFC- 















13^6 1318 HOB Hia . 15C8 IS18 

TIME PERIOD 



IBdS IB'18 



Figure 46. 



Same as Figure 37 except for Vertical Vorticity 
Advection . 
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Same as Figure 37 except for the Tilting Term 
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8. Same as Figure 37 except for Frictional 
Dissipation. 
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Figure 49. 



Same as Figure 37 exceot 
Positive Values Indicate 
Negative Values Indicate 



for Budget Residuals. 
Vorticity Excess. 
Vorticity Deficits. 
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Figure 50. FGGB 6-h Radius 6 Vorticity Budget Fields for 
(A) Time Tendencies (B) Buaget Residuals. 
Period 1203 Refers to 0000-0600 GMT 12 January 
and is Valid for 0300 GMT. Contours as in 
Figure 37. 
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Figure 50. FGGB 6-h Radius 6 Vorticity Budget Fields for 
(A) Time Tendencies (B) Buaget Residuals. 
Period 1203 Refers tc 0000-0600 GMT 12 January 
and is Valid for 0300 GMT. Contours as in 
Figure 37. 
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Figure 51. FGGE 6-h Stability Trace. Units and Times as in 
Figure 28. 
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